Integrability and Ergodicity of Classical Billiards in 

a Magnetic Field 

N. Berglund, H. Kunz 
Institut de Physique Theorique 
Ecole Polytechnique Federale de Lausanne 
PHB-Ecublens, CH-1015 Lausanne, Switzerland 
e-mail: berglund@iptsg. epf 1 . ch, kunzOeldp . epf 1 . ch 

Revised version May 1995* 



Abstract 

We consider classical billiards in plane, connected, but not necessarily 
bounded domains. The charged billiard ball is immersed in a homogeneous, 
stationary magnetic field perpendicular to the plane. 

The part of dynamics which is not trivially integrable can be described 
by a "bouncing map". We compute a general expression for the Jacobian 
matrix of this map, which allows to determine stability and bifurcation values 
of specific periodic orbits. In some cases, the bouncing map is a twist map 
and admits a generating function. We give a general form for this function 
which is useful to do perturbative calculations and to classify periodic orbits. 

We prove that billiards in convex domains with sufficiently smooth bound- 
aries possess invariant tori corresponding to skipping trajectories. Moreover, 
in strong field we construct adiabatic invariants over exponentially large times. 
To some extent, these results remain true for a class of non-convex billiards. 

On the other hand, we present evidence that the billiard in a square is 
ergodic for some large enough values of the magnetic field. A numerical study 
reveals that the scattering on two circles is essentially chaotic. 

Key words: billiards, magnetic field, twist map, integrability, 
adiabatic invariant, ergodicity 

1 Introduction 

We consider the classical motion of a particle of mass m and charge g in a plane 
domain Q. A homogeneous, stationary magnetic field B perpendicular to the plane 
makes the particle move on arcs of Larmor radius ^. Whenever it encounters the 
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boundary, the particle is reflected specularly. This problem was first considered by 
Robnik and Berry [|R^). 

One point of interest in such models is the problem of integrable versus ergodic 
behaviour of Hamiltonian systems. In zero field, we know examples of integrable 
billiards (elliptic or rectangular boundary) as well as ergodic ones (dispersive Sinai 
billiards, Bunimovich stadium, see for instance ||KT|| ). Since it is unlikely that 



integrability is stable with respect to perturbation by a magnetic field, i? is a natural 
parameter for studying the transition from order to chaos. One question is whether 
the billiard can become globally ergodic in strong enough field. 

Another motivation for studying magnetic billiards is connected to the problem of 
quantum chaos. Recently, it has become possible to create mesoscopic systems where 
the electrons' motion is essentially ballistic [[LRPW|| . Behaviour of some macroscopic 
observables, like the susceptibility, can be surprisingly complicated ||RKC|| , and may 
be related to classical dynamics. In order to be able to apply semi-classical methods 
like trace formulae, it is desirable to have a good knowledge of classical periodic 
orbits and their stability, at least perturbatively. 

This paper is organized as follows. In section 0, we give a more precise definition 
of the billiard flow and boundaries considered. The interesting part of dynamics 
can be described by a bouncing map T, defined in section ^, where we also give an 
exact expression for the Jacobian matrix of T. In some cases, T is a twist map and 
admits a generating function G. In section ^, we give an exact expression of G, and 
we discuss its physical interpretation and its practical applications. 

Our main analytical results are contained in section ^, where we consider bil- 
liards with smooth, convex boundaries. Using perturbative techniques, we analyze 
different quasi-integrable limits. In section p.2\ , we prove the existence of invariant 
curves for all values of the magnetic field. These quasi-periodic trajectories corre- 
spond physically to diamagnetic currents along the boundary. Our KAM type result 
generalizes the well known theorem of Lazutkin on the existence of caustics near 
the boundary, in the zero field case. The magnetic field however breaks the symme- 
try between the forward and backward skipping orbits, and at intermediate values 
of the field, only one kind of skipping orbit is present. In section we analyze 
the strong magnetic field limit. We compute an expansion of the bouncing map in 
powers of the Larmor radius and construct an adiabatic invariant on a time scale 
growing exponentially with the magnetic field when the boundary is analytic. At 
low order, this adiabatic invariant coincides with the one derived by Robnik and 
Berry. In order to construct the invariant, we derive a theorem on a class of maps 
of the annulus, which may be of interest in a broader context. Section ^ is devoted 
to the proof of it. 

In section ^, we study billiards which may show an important chaotic compo- 
nent, because of singularities in the boundary (polygons) or because of a negative 
curvature. It turns out that the properties of billiards in convex domains can be 
to some extent generalized to billiards with concave boundary. In general it seems 
that billiards in a magnetic field are of the mixed type (invariant tori and chaotic 
components coexist). A possible exception is the case of the square at some par- 
ticular values of the magnetic field, where we present an analytical argument and 
numerical evidence for a completely chaotic motion. 



Figure 1: Examples of billiard domains Q, with tangent and normal vectors to the 
boundary dQ: (a) bounded domain, (b) unbounded domain. 



2 Definition of the billiards 

Let Q be a connected domain in IR^. We assume for simplicity that the boundary dQ 
consists of simple, closed, piecewise curves of total length \dQ\, although some 
of our results may be extended to more general cases. Several of them however need 
a higher degree of differentiability, which will be clearly indicated on the spot. 
The boundary is parameterized with the curvilinear abscissa or arclength s: 

x(s) = (X(s),F(s)), ds^ = dX^ + dY\ se[0,\dQ\). (1) 

The unit tangent and normal vectors and the (signed) curvature are given by 

t(s) = = (cosr(s),sinr(s)), 

n(.) = {-Y'{s),X'{s)), 

dr 1 
k{s) = - = X'{s)Y"{s)-X"{s)Y'{s) = -^^. (2) 

Parameterization is chosen in such a way that n is always oriented towards the 
interior of Q (see fig.[T|). In this way, the curvature is positive for a convex boundary. 
We suppose k(s) to be defined and continuous everywhere but on a set E2 of punctual 
values of s. The vectors t and n are undefined on a set Ei <Z E2. 
Inside Q, the billiard flow is given by the Lagrangian 

L(x, x) = imx^ + q (i| A(x)) , A(x) = {-\yB, \xB). (3) 

The resulting motion is simply circular uniform, with Larmor radius 

mv y/2mE , > 



where speed f = |x| and energy E are constants of motion. We adopt the sign 
convention qB < 0, which imphes that the particle turns counterclockwise, and 
we take some characteristic length of the billiard as a length unit, so that n can 
be considered as a dimensionless parameter. Thus, increasing B is equivalent to 
decreasing E or increasing the size of the billiard (which can be viewed as taking 
the thermodynamic limit). 

Each time the particle hits the boundary, the velocity changes according to the 
law of specular reflection v' = v — 2 (v|n) n, which is well defined for almost every 
point of the boundary^]. 

The phase space M = Q x can be divided into two disjoint sets Mi and M2. 
Ml consists of all the orbits that never touch the boundary. It corresponds to an 
integrable component of the motion and may be empty. In non-zero magnetic field 
{fi < 00), the orbits of M2 hit the boundary an infinite number of times. It is thus 
natural to study the dynamics in M2 by means of a bouncing map. 



3 The bouncing map and its Jacobian matrix 

Let us consider the trajectory between two successive collisions with dQ, occurring 
at Pq = (-^(■So)) ^(■^o)) and Pi = (X(si), Y{si)). The trajectory is an arc 7 of center 
O, radius /i and angle ip (see fig.@). 

We call 6i the angle between the arc and the boundary at Pj, and Ui = — cos 6i, 
i = 0, 1. Quantities s and u are the Birkhoff variables and the bouncing map T is 
defined as 

T : (so,Uo) ^ (si,Mi). (5) 

If i is the length of the chord PqPi, and x is the angle between the chord and 
the arc 7, simple geometry shows that 

^ = 2x, smx = — . (6) 

In general, there may be two trajectories with supplementary x for a given £ (see 
also fig.|TB|). This is a characteristic magnetic field effect. 

Generalizing Birkhoff 's technique (Q, p. 173, see also |[KT|| ), we find 

Proposition 1 If Sq, Si ^ E2, the Jacobian matrix DT ofT is given by 

dsi Koi cos X — sin(6'o + 2%) 

dso sin 9i 

dsi £ cos X 

duo sin 9 1 sin 6*0 

dui sin(^o + 2x) sin(6'i + 2%) — sin Oq sin 61 

dso i cos X 

—kq sin(6'i + 2x) — fi^i sin(6'o + 2%) + kqHiE cos x 

dui Kid cos X ~ sin(6'i + 2%) 



dun sin 6'i 



(7) 







"'^ Subtleties may occur in the corners or if (v|n) = and p < at the point of coUision, but 
they are of Uttle importance for the foUowing. 



o 



Figure 2: The trajectory between two successive bounces is an arc 7 of radius /x, 
angle if) and extremities Pq and Pi. The Jacobian matrix of the bouncing map T 
depends on the length ^ of the chord, the angle x between chord and arc, the angles 
6*0, 6*1 between arc and boundary and the curvature at Pq, Pi- 



where Hi = K,{si). 

For completeness and consistency of notation, we give a proof in appendix 
We noticed that an alternative proof of this formula was provided by [|MBGt|] . All 
the quantities appearing in (|^ can be easily expressed as functions of sq and si (see 
([76|) and ([77|) ). It is in general much more difficult to solve the relation mo(so, si) 
with respect to Si, in order to obtain DT as a function of zq = (so,Mo)- 

A first observation is that DT has unit determinant, and hence the Birkhoff 
variables are conjugate, although momentum and velocity are not collinear in a 
magnetic field (see section 

A second observation concerns the element We see that its sign depends 
only on x- Suppose that the shape of Q is such that it cannot contain any arc 
of radius fj. and angle larger than tt. Then we have ^ > 0, and uq is uniquely 
defined for given sq and si. In such a case, T is called a (symplectic) twist map. 



Twist maps, especially when continuous, have a lot of properties (see e.g. [[Me|| for 
a review), some of which we will discuss in the next section. 

The Jacobian matrix allows us to compute Liapunov exponents and to determine 
stability of periodic orbits. If we define for z = (s, u) the stability matrix 

Sn{z) = DT{T''-h) ■ ■ ■ DT{Tz)DT{z), (8) 

then the Liapunov exponents are given by 

A±(z) = ±lim^lnTTS:{z)Sn{z). (9) 



If z belongs to an orbit of period n, the orbit is hyperbolic and unstable if 
|Tr5'„(2;)| > 2, parabolic if |TrS'„(z)| = 2 and elliptic if |TrS'„(2;)| < 2. In the latter 
case, the orbit is stable unless resonance occurs. In section ^, we will give examples 
of how to apply the formula (0) to analyze stability of specific families of periodic 
orbits. 



4 Generating functions 

Here we call generating junction of the continuous map T a function G(so, s\) 
such that 

dG = uqcIso — Uidsi. (10) 
It is easy to check the following properties ||Me|| : 



1. If T is an area preserving twist map, it admits a generating function, unique 
up to an additive constant, given by 



Uoi^,v)d^ - M^,v)dv- (11) 



2. If G is C^, the map T generated by G is always area preserving. It is a twist 
map iidl,^G>Of\. 

3. If G is C^, M is a constant of motion iff G{so, Si) = g{si — sq). 

In zero field, G is known to be the length of the chord i{so,si). For magnetic 
billiards, we found that G depends also on an area associated to the trajectory, a 
feature appearing apparently in all problems involving a magnetic field: 

Proposition 2 Suppose that Q is bounded and that T is a twist map. Then, for 
So,Si ^ E2, the generating function is given by 

G = C + -S, (12) 
/i 

where C is the length of the arc and S is the area between the arc and the boundary 
(see fig^. 

We give the proof in appendix ^. For practical purposes, it is useful to write G 
in the form 

G(so, si) = -A{so, s,) + B^{i{so, si)), (13) 

where A, the area between the chord and the boundary, does not depend on the 
magnetic field, and 



B,{e) = /iArcsin ( A ) + ^ /i _ (14) 



2fil 2 V 4/x 



can be degenerate if dl^^^G = 0. 



Figure 3: When T is a twist map, its generating function G can be expressed (a) 
as a function of the length £ of 7 and the area S between 7 and dQ, or (b) as a 
function of the length i of Pq-Pi and the area A between Pq-Pi and dQ. 



is equal to C minus the area between the chord and the arc, divided by /i. 
As an example, let us consider an elliptic boundary, of equation 

ds 



X = (Acosy^, siny^), —— = C{(p) = \/ cos^ (p + X"^ sin^ cp. (15) 

dip 

If the magnetic field is low enough (in fact, if yU > A^, as we shall see in section |^), 
T is a continuous twist map, and 

A 



G = 



1 , 
- - sm{2(p. 



+ P^(2sin¥;_C(^+)), (16) 



where ip± = ^{ipi±(pQ). Note that for a circular boundary (A = 1), C = 1 and G de- 
pends only on ipi — ipQ. Thus, the billiard becomes integrable, which is geometrically 
obvious. 

In order to give a physical interpretation to G, let us recall that the momentum, 
which is canonically conjugate to the position, is given by 

dL X 1 

p = — — = mx + gA = (nix qBy, my H — qBx). (17) 

ax 2 2 

Thus, it would have been physically more natural to use, instead of the tangent 
velocity m, the tangent momentum 

P=-— (P|t(.)) =u + ^{X{s)Y\s) - X'{s)Y{s)). (18) 

Indeed, p and s are conjugate since they admit as a generating function the reduced 
action along the arc 7: 

F = I (p|iix) = mv f ds + — {ydx — xdy) =^ dF = mv{podso — pidsi). (19) 

It is however more convenient to use u instead of p. Green's theorem implies that 
the generating functions are related by 

mvG = F + -qB f ydx — xdy, (20) 

2 J'JQ 



where 7q is the piece of boundary connecting Pq to Pi. 

One apphcation of generating functions is in perturbation theory. Suppose the 
map depends on a parameter e (controlhng the magnetic field or the shape of the 
boundary), such that the behaviour of Tq is known (e.g. integrable). Approximating 
Tg by its expansion TQ + edeTe\£=Q does in general not lead to an area preserving map, 
whereas the map generated by G{sq, Si, Q)+edeG{sQ, Si, £:)|£=o is always conservative. 

For example, the generating function for an elliptic boundary close to a circle, 
given by (|T5|) with A = 1 + e, is 



G 



Go{^- 
G'i(y?+,y?- 



1 

^^ 
1 
/i 



V- 



V- 



+ 2 sin(2(^_ 



+ 2 sin(2(^_ 



+ P^(2 siny9_ 



+ 2 sin sin ^p^ 



sin^ (p. 



(21) 



Another use of generating functions is in searching periodic orbits. If we define 
the n-point generating function 



G^'^\so, si, . . . , s„_i) = G{so, si) + G{si, S2) H h G{sn^i, sq) 



(22) 



then the law of specular reflection implies that every periodic orbit containing no 
points with s G £'2 is a solution of 



dsr 



dsi 



ds. 



0, 



(23) 



n-l 



which is a system of n nonlinear algebraic equations of the n variables sq, . . . , Sn-i- 
It is convenient to "lift" the periodic variable s to the real line. An (m, n) 
periodic orbit is then defined by Sn = sq + m\dQ\ ,Un = uq, and its frequency 
u = lim,fc_,oo ^ is equal to ^. The possible frequencies belong to an interval I^, 
depending on the behaviour of the boundaries u = ±1 of the phase cylinder. First 
results on existence of periodic orbits for continuous area-preserving twist maps are 
due to Poincare and Birkhoff. Powerful developments were achieved by Aubry and 



Le Daeron, Mather, MacKay and Meiss, and Katok (see ||Me|| and references therein 
for more details): 



1. For every m, n, — G /^j, there is at least one 



periodic orbit which is 



"maximizing". This means that every finite orbit segment (s^, . . . , s;), / > 
/c + 2 is a global maximum of J2''j^kG{sj, Sj^i) with respect to variations of 
Sfc+i, . . . , S;_i. In particular, (sq, . . . , s„_i) is a global maximu mQ of If 
the maximum is nondegenerate, the orbit is hyperbolic. 

2. For every m, ra, ^ G 1^^, there is at least one {m,n) periodic orbit which is 
"maximin". This means that the Hessian matrix of I]j=fc ^^(■Sj, sj+i) has one 
single positive eigenvalue. The orbit is either elliptic or inverse hyperbolic 
(TrS-, < -2). 



^Second variations of G can be computed using (|o|) and (pl|). 



3. Every orbit on a rotational invariant circle is maximizing. For every irrational 
uo G luji there is a maximizing quasiperiodic orbit of frequency uj. Its closure 
is either an invariant circle, or an invariant Cantor set. This result is in some 
sense stronger than KAM theory, since it shows the existence of quasiperiodic 
orbits for twist maps that are not necessarily nearly integrable. 



To summarize: if the bouncing map is a continuous twist map, its generating 
function provides a useful tool to do perturbative as well as variational calculations. 
In particular, we have a lower bound on the number of periodic orbits of period n 
(namely twice as many as there are integers m coprime with n such that — G J^^), 
whose stabihty can be related to the second derivative of the generating function. 

As we shall see in the next section, an important class of maps satisfying the twist 
property is given by bouncing maps of billiards with smooth, convex boundaries in 
low magnetic field. Many billiards' bouncing maps however are either not twist or 
discontinuous. Nevertheless, generalized versions of the generating function can be 



useful tools in some of these cases as well, as we will see in sections 5.3 and 3.2 



5 Smooth convex boundaries 



5.1 Qualitative behaviour 

We say that dQ is smooth and convex^ if p{s) is a smooth function (in a sense that 
we will have to specify) and is bounded by positive constants: < pmin < p(s) < 
Pmax < 'is. For example, for the ellipse (p!5D, pmin = and pmax = A^. 

Robnik and Berry ||RB|| proposed to classify the dynamics by comparing to 
Pmin and Pmax- If wc cousidcr specifically the twist property, by looking at trajecto- 
ries with fixed Sq and varying uq (fig.§ and |), we obtain the following "curvature 
regimes" : 



1. If p > Pmax, the function si(so,mo) is increasing in mq, so that T is a twist 
map. The curve {T(so, Mq), —'\-<uq< 1} turns exactly once around the phase 
cylinder, with lim^.o^±i T(so, mq) = {sq.Uq). Thus, = [0,1]. 

2. If pmin < p < Pmax, discontinuities in the map occur when the trajectory 
becomes tangent to the boundary. We have lim„(,^+i T(so, Mq) = {sq.Uq), but 
this is not necessarily true for uq ^ —1. In appendix 0, we illustrate the 
construction of the discontinuity lines of the map. 

3. If p < Pmin, the function Si{sq, Uq) is first decreasing in uq, reaches its minimum 
when X = ^/2, and increases again. We again have lim„Q^-|-i T(so, Mq) = 
(scMo), so that there are exactly two trajectories with supplementary x ^oi 
given So, si. The map is no longer twist. One can expect that the frequencies 
belong to an interval = [a;min, 0], where tUmin ~ p- 



"'For brevity, we omit the word "strictly". 



Figure 4: Trajectories in an elliptic boundary, with constant Sq and varying uq: 
(a) fi > pmax, (b) pmin < ^ < p^ax, (c) p < Pmin- Trajectories with Mo near -1: 
(d) trajectory with caustic for p > Pmax, (e) chaotic trajectory when pmin < p < 
Pmax, (f) trajectory with caustic for p < Pmin- 



Figures ^ and ^ illustrate this behaviour in the case of an elliptic boundary^. 
We see that there always exist invariant curves near u = +1, corresponding to 
trajectories with caustics (the "whispering gallery modes"). For p > p^ax and 
P < Pmin, such curves also exist near u = —1. However, when pmin < p < Pmax, 
the tori near u = —1 are replaced by a chaotic region. This behaviour can be 
heuristically understood by noting the following two points: first, discontinuities of 
the map due to tangencies may destroy invariant curves in an analog manner as do 
the discontinuities due to corners of the boundary [ DR |; second, the Jacobian matrix 
(0) diverges like + causing strong dispersion of nearly tangent orbits, which 

can be a source of positive Liapunov exponents. 

Robnik and Berry have computed an adiabatic invariant, which tends to be 



'''The billiard in an ellipse has the particularity to possess two integrable limits. In the circle 
limit A = 1, u is a constant of motion, as we noted in the previous section. In the zero field limit 
/i = oo, the product of the angular momenta with respect to the foci is a constant, which can be 
written L = C'^{ip){u'^ - 1) + 1. 




Figure 5: Structure of phase space in the {ip, u) plane for the eUipse. Image of the 
hue ((/? = (/?*,-! <u< 1) : (a) /i > pmax, (b) Pmin < /i < Pmax, (c) fJ, < Pmin- Typical 

phase portraits: (d) > Pmax, (e) Pmin < A* < Pmax, (f) A* < Pmin- 



conserved if 



l + KflU 



< 1: 



K{s,u) 



(1 - ^2)3/2 



(24) 



and gives a very good approximation of the invariant curves for m ^ ±1. As it 
turns out, in the strong field hmit fin 1, K tends to be conserved for all orbits, 
so that this hmit seems to be integrable. In the foUowing, our aim is to dehneate 
the vahdity of these statements and make them more precise by using different 
perturbative approaches. 



5.2 Near the boundary 

We recall that for Euclidean billiards, the existence of invariant curves and caustics 
near the boundary was first demonstrated by Lazutkin [Q], who had to assume 
a high degree of differentiability. Douady |po|| reduced the required degree to 6. 
On the other hand, Hubacher H] proved that no caustics exist near a boundary 



whose curvature is discontinuous, and Mather ||Ma|| found that if the curvature of 



the boundary vanishes at one point, then invariant curves are totally absent. 

We are now going to investigate the behaviour of the bouncing map of billiards 
in non-zero magnetic field near their boundary, i.e. for small sin 6*9. In this way, we 
will be able to apply KAM theorems to show the existence of invariant curves. The 
results are summarized in theorem 1 at the end of this section. 

If we write z{s) = X{s) + iY{s), and Zi = z{si), Ti = r(si), we have 



''TO 



zi-zo= I e'^'^'^ds= / p(r)e'"rfr. (25) 



Since the Zi also belong to the arc of trajectory, which has tangent directions 
To + 9q at Zq and ti — 61 at zi (see figj^), we have 



_ ^^ 

Zl — Zq — — 
1 



gi{ri-ei) _ gi{ro+eo) 



(26) 



Introducing A = e '^°(-2i — 2^0) ^'Hd 5 = ri — Tq — ^1 + ^o; this can be rewritten 

A - /ie-^^°^^ = -2/isin^o, (27) 
i 

which is equivalent to the system 

Re A ^ sm6 , . ^ \ 

— /icosyo — + /i(coso — 1) = — 2/i 

sin 9o sin Uq 

ImA sin 6 cos 6 — 1 

— 2^ +/^— ^ + /^cos^o— = 0- (28) 
sm uq sm uq sm uq 

If the boundary is C'', this is a system of C^~^ equations in the variables So,6q,si, 61 
that we would like to solve with respect to Si, ^1. 



Writing si = sq + sin 9ofi<J and 9i = Oq + sin 6*0/17, we obtain for A; > 3 



A = sin 6ojj,a + - sin^ O^^'^a'^ k{sq) + 0(/i^ sin'^ 6*0) 



.2_2. 

.2 „• 2 , 



5 = sin^o/u[«:(so)(^-7] +0(/i^sin^^o)- (29) 
Replacing this in (EBf) finally leads to the system 



cr — /icos^^o ['^(so)c'" — 7] + 0(sin6'o) = —2 
(J^k(so) + 2 [K(so)cr - 7] - /icos6'o [k(so)o- - 7]^ + O(sin^o) = 0, (30) 



which has the solution 



_2 

a{so,9o,fi) = — — + 0(sin6lo) 

1 — /i COSUoK[So) 

7(so,6'o,/i) = 0(sin6'o). (31) 



Since the Jacobian of ( pOD evaluated at this solution is —2 + O(sin^o)) the implicit 
function theorem implies that for small sin6'o, the bouncing map is C*^"^ and has 
the form 

2/isin6'o / . n N / 1 lo^iN 

si = So-- — -— + o(sm6'o) (mod |dQ|) 

1 — ncostloi^{so) 

61 = eo + o{sm9o). (32) 

However, this formula is only valid if we are able to check two properties: firstly, 
the approximation must be well defined, i.e. the denominator in the first equation 
should not vanish. Secondly, the solution ( pT] ) has to be the "right" one, that is the 
first intersection of the trajectory with the boundary, following the motion of the 
particle. We have to distinguish between the following cases: 

1. Near u = +1: Writing 6 = n — t], (B^I) becomes 



si = So- r]o + o{r]o) (mod |aQ|) 

1 + Ij,k{so) 

Vi = Vo + oivo)- (33) 

Here the denominator can never vanish. Moreover, this solution corresponds 
geometrically to a short skip of the particle backwards along the boundary, so 
that it certainly describes the first encounter of the trajectory with dQ. 

Note that this equation is well defined even when the curvature is allowed to 
vanish. Only when /i = 00, the first equation becomes si = so — -;^i^rio + o{rio), 
and diverges when k — 0, which is coherent with Mather's result. 

The change of variables (p = s + /ur(s), r = 2^1] is well defined and transforms 
the map into 

ipi = (po-ro + o{ro) (mod \dQ\ + 27r/i) 

n = ro + o(ro). (34) 



2. Near u = —1: The map has the form 

si = So- ^ , ^0 + 0(^0) (mod |9g|) 

Oi = Oo + oiOo). (35) 

This time, we have to be careful with the denominator, which may vanish. 
The approximation ^ -C 1 is only valid in two cases: 

(a) If /i > Pmax(l + e), e > 0, then in the variables ip = /ir(s) — s, r = 2fi9, 
the map becomes 

V^i = V^o + ^0 + 0(^0) {mod 27T jj, - \dQ\) 

n = ro + o(ro). (36) 

We again got the right intersection, because the particle is skipping for- 
ward along the boundary (fig.^). 

(b) If /i < Pmm(l — £), then in the variables (p = s — ^t{s), r = 2fi9, the map 
reads 

<fi = (po-ro + o(ro) (mod \dQ\ - 2nfi) 

ri = ro + o(ro). (37) 

This time, the particle is starting with a forward glancing velocity, but 
reaches the boundary behind its starting point (fig.^). It is not obvious 
that there is not another intersection of the trajectory and the boundary 
in between. Luckily, lemma 3 in appendix |^ shows that this cannot 
happen, because any circle of radius < cuts dQ at two points at 
most. 



The maps (^4]), (^) and (|37|) can now be treated by KAM-type theorems ( ||Mo|| , 
p. 52 or [ po| , p.III-8). In this way, we obtain the following result: 

Theorem 1 Consider a magnetic billiard in a domain with boundary, k > 6, 
such that < pmin < p(s) < pmax < CO. Lct 6 G JR\ and assume that one of the 
following is true 

1. 0<p,<oo,r] = n-6,L = \dQ\ + 2np,, a = -1; 

2. Pmax(l + e)<pi<oo,ri = 9,L = 2tx^ - \dQ\, a = +1; 

3. < p < p^i„(l -e),7] = e, L= \dQ\ - 27rp, a = -1. 

Then, there exists > (depending on fi and k) such that, if uj [0,£:i) and 
satisfies the diophantine conditions 



U! p 

L~q 



> 7g-^ (38) 



for some 7, G H!^ and for any - G (Q, there is an invariant curve of the form 



V = ^ + ^(0, (39) 

where U,V e C\ V{i + L) = 1/(0 + \dQ\ - L, + L) = U{0- The induced map 
on this curve has the form 

i^i + cru. (40) 



This theorem confirms our observation that invariant curves exist in three cases: 

1. Near m = +1, {9 = -k), for all values of the magnetic field (see the upper parts 
of fig.^-f). They correspond to backward skipping trajectories, the curvature 
of which is opposite to the curvature of the boundary. 

2. Near m = — 1, (0 = 0), in weak magnetic field (/i > Pmax(l + ^)-, see the lower 
part of fig.pi). They correspond to forward skipping trajectories (fig.§i), 
which are curved away from the boundary. 

3. Near u = — 1, (6* = 0), in strong magnetic field (/i < Pmin(l — see the lower 
part of fig.^f). They correspond to backward skipping trajectories (fig.|f). 

For intermediate values of the magnetic field, there seem to be no tori near u = 
— 1. Although we have a qualitative understanding of the origin of this phenomenon, 
we are not able to prove the existence of a stochastic component of positive measure 
in this region. This remains an open problem. 

On the other hand, our theorem implies that magnetic billiards in smooth (i.e. 
C^) convex domains can never be ergodic. Unlike in zero field, this remains true 
even when the curvature of the boundary is allowed to vanish. In fact, equation 
(P3|) suggests that invariant curves exist near u = +1 even when the boundary has 
slightly concave parts, i.e. k{s) > —Kq, Kq > 0, and for fi < — e. In that case, 
however, the tori near m = — 1 do not necessarily survive. We will come back to that 



point in section 6.2. 



5.3 The strong field limit 

The limit /i ^ is in some sense singular: indeed, for small /i, the distance between 
successive collisions with the boundary is of order /i, so that one could get the 
impression that the particle follows the boundary more and more slowly when /i 
decreases, and sticks to the wall when /i = 0, violating the conservation of energy. 
Of course, this is only an artificial effect, due to the fact that we use the bouncing 
map instead of the flow to describe the dynamics. Indeed, the "time of flight" 
between successive collisions is also of order fi (more precisely, it is equal to fiip, 
where ip is given in fig.0), so that during a given time interval, the particle hits the 
boundary 0{l/fi) times, traveling a finite distance. 

Thus, we have to be careful when we apply perturbation theory. It is possible to 
use (|2^) to analyze the limit /i — > 0. However, we give in appendix an alternative 



approach that we find more instructive. The idea is to show that although the map 
is not twist, it can be described by generating functions, in fact two of them, with 
suitable matching at x = 7r/2. These functions are of the form G^{so, si, fi) = 
±figo{a) + /i^f^f (cr, So, yu), where a = Replacing this in ([T0|) , we see that a 

and u are "slow" variables which evolve on a much longer time scale than s. This 
expresses the geometrical idea that for a short skip, the boundary is close to an arc 
of circle, for which a and u are constants. It is however necessary to use the variable 
6 instead of u if we wish that the map be smooth at the boundaries of the phase 
cylinder. Finally, we obtain the following result: 

Proposition 3 Consider a magnetic billiard with a convex, boundary, k > 3. 
For small enough fi, the bouncing map is C^~^ and has the form: 



$1 = So — 2/i sin 6*0 + /i^ sin 6*0 a(so, 6*0, yu) (mod |<9(5|), 
Oi = 9o + fi^sm'^9ob{so,9o,n). 



(41) 



The functions a E and b E ^ are uniformly bounded for sG IR, 0<^^<7r, 

\dQ\-periodic in s, and admit the expansions 



a(s, 6, /i) 

6(s,e,/i) 
where the first terms are 

ai 

ho 
bi 



k—3 

J2a,{s,e)fi' + 0{fi' „ 

i=0 
3=0 



.k-2\ 



,k-3\ 



(42) 



—2 cos9k{s), 

[(1 + 2 cos 2e)K{sf - sin 29 k (s) 

3 ^ 

4 2 
- cos 9 k,{s)k,'{s) — -sin6'fi;"(s). 



(43) 



Equation (^T]) shows that the bouncing map in strong magnetic field behaves 
like a perturbed integrable map. To zeroth order in /x, it reduces to identity. To 
first order in /i, it still has the integrable form 9i = 9q, si = So + ^Q{9q), where 
Q{9) = —2 sin 9. However, in contrast with usual integrable systems, the frequency 
Q is multiplied with the small parameter fi. As for the factors sin^ occurring in 
(PI), they assure that the boundaries 9 = 0, n are fixed. 

The behaviour when — > can be understood in the following way: fix some 
positive /iQ and define s = ^Qip, d{ip,9,fi) = a{fioip,9, fi), b{(p,9,fi) = b{fio{p,9, fi). 
Consider the 2-parameter family of maps T^q.m- 

<fi = V'o - 2sin6'o + /isin6'oa(v9o,6'o,/i) (mod ■^^^), 

9^ = ^o + /i'sin2^o%o,^o,/i)- (44) 



The map T^o,o is integrable, 9 being a constant of motion, and it takes about 
n{9) = 2j(f^ig iterations for Lp to turn once around phase space. T^o,mo i^ equivalent 
to (|4l|) , and one can expect that when fiQ is small, some features of the integrable 
map remain, in particular the particle completes one turn after roughly n{d) bounces. 
Using (|77D and (0), we obtain that ip = 2x = — 29q + 0{fi), so that the time 
necessary for one revolution is of order t ~ n{9)noip ~ \dQ\ which does not 
depend on /i to lowest order. Note an important difference between the skipping 
regimes m ~ +1 and m ~ — 1: in the first case, t ~ \dQ\, whereas in the second 
case t ~ \dQ\ | diverges when m — >■ — 1. The variation of 9 during t is of order 
n{9){9i — 9q) ~ /iosin^, which shows that our approximation is consistent, since 9 
is an adiabatic invariant on the time scale t. 

We have the choice between two different perturbative techniques to make these 
observations mathematically precise. If we exclude some neighborhood of ^ = |, 
where the frequency VL{9) = —2 sin 9 has a vanishing derivative, we can apply Moser's 
theorem ( | [Mo|| , p. 52): if the boundary is C^, k > Q, then there exist invariant curves 
for every frequency satisfying some diophantine condition^. In this way however, we 
cannot describe trajectories starting with a very small tangential velocity {9 ~ |). 
In order to do this, we will apply another technique to study (^T]), coming from 
adiabatic theory. We will prove in section [7| the following theorem on a class of 
maps of the annulus: 

Theorem 2 (An adiabatic theorem for maps) 

Let the map T : {ipo, Iq, /i) i-^ (ipi, Ii), 



h = Jo + /i^c(/o)^/3(v?o,^o,/^), 



(45) 



be defined on the set E = {{{p,I,fi) Iv^GlR, 0</<l,0</i< fi*}. Assume that 
a and (3 are periodic functions of ^9 with period 1, c(0) = c(l) = 0, and 0) = 

Then, there exists a change of variables ^ {ip = <f + ^f{(f,I,fi),J = 

I + fic{I)g{ip, I, jj)), preserving the square [0, 1] x [0, 1], where f and g are periodic 
with period 1 in (f, such that 



1. IfTis analytic in a complex neighborhood of E, then 

^1 = V'o + M<^o) [^(^o,yu) +e"^/^''a(^/'o, Jo,/u) 
Jl = Jo + /iM'/o)'[e(Jo,/i)+e-i/^'^^(V'o,Jo,/i) 



(46) 



2. If moreover T preserves the measure c{I)p{ip, I , fi)d(pdl , where p{ip,I,0) = 1 
then 

Jl = Jo + e-^/^M'/o)'^(^o, Jo,/"). 



(47) 



3. If a, (3 and Vt are C'^ and c is C^^"" in E, then 

Jl = Jo + fi^c{Jof Q{Jo,fi) + fi''/3{iJo,Jo,fi) 



(4J 



'One advantage of Moser's theorem is that the map need not preserve an annulus. 



respectively 

Ji = Jo + fx'+^ciJoyPii^o, Jo, /i) (49) 
ifT preserves the measure c{I)p{ip, I , fi)dipdl . 



Equation (|45|) describes a map from the annulus < / < 1 into itself (the 
choice of 0, 1, is arbitrary and can be changed by scahng). Indeed, the factor c(/), 
which vanishes at / = 0, 1, assures that the boundaries of the annulus are fixed, 
and for small enough /i, the interior of the annulus is invariant. In (|45|) , the terms 
containing the phase cp, which make the map non-integrable, are of order /x^. The 
theorem states that a suitable change of variables decreases this order to fi^'^'^ if the 
map is C^, and to the exponentially small order e~^^^^ if the map is analytic. 

The method of the proof is similar to that used by Nehoroshev and Neishtadt 
|A2|| , p. 163 and references therein) in the case of differential equations. The 



see 



basic idea is to make successive changes of variables that each decrease by one unit 
the order of the terms containing the phase ip. If the map is C'^, one can repeat 
this procedure k times. If it is analytic, divergence prevents us from applying the 
procedure infinitely often, but it is possible to do it a large number of times, of order 
^, so that the terms containing the phase become exponentially small. 

We can now apply the theorem to the bouncing map ( ^T] ) of billiards in strong 
magnetic field, where c{6) = sinO. Because of proposition 1, the map preserves the 
measure sin 6dsd6, and can thus be transformed into the form (^) or (^91). If we 
write Ji = J{si,9i), we obtain 

Corollary Consider a magnetic billiard in a convex domain. If the boundary is 
C^,k > 3, then there exists a function J^'^^{s,6) such that Ji'^ = Jq'^ + 0{ix^~^). 
If the boundary is analytic (in the sense that n{s) can be analytically continued to 
a complex neighborhood of the real axis), then there exists a function J{s,9) such 
that = Jo + 0(e-i/^'')- In any case, this function can be written in the form 
J = 9 + fisin 9g{s, 9, /i) . 

If the boundary is C'^, A; > 3, this result implies that after n bounces, we have 
Jn = Jo + nO{fi^~^), and hence if n = 0(/i^~^), J„ = Jo + 0(/i). In other words, the 
quantity J is a "quasi-invariant" or "adiabatic invariant" which varies of an amount 
of order /i for 0(/i^~^) bounces, i.e. during a time of order In the analytic case, 

J varies on a time scale growing exponentially with the magnetic field. In contrast 
with KAM theorems, which show the existence of exact invariants for some initial 
conditions, the adiabatic theorem shows the existence of approximate invariants for 
all initial conditions. 

In fact, the theorem not only shows the existence of quasi-invariants, it also 
enables us to compute them up to the first orders. For example, after 2 changes of 
variables (given by (|5^), (|53|) and (^ID), we obtain 



J^^\s,9,fi) = 9 + fism9 



12 

-k(s) + -fi cos 9k{s) 



(50) 



To zeroth order in /i, we recover that J*^'^^ = ^ is an adiabatic invariant on a time 



2iu ^0 



2ir 



Figure 6: Adiabatic invariant in the strong field limit: (a) some orbits of the billiard 
in an ellipse for = 0.3, A = 1.5 (i.e. pmin = |), (b) level lines of the adiabatic 
invariant (0), which are given by 6* = J — /i sin J + |/i cos Jk(s)^ + 0(/x^). 



scale of order 1. To first order in /x, J*^^^ corresponds to the invariant p^ ) of |P^B|| p], 
which is valid for a time of order As shown in figure ^ the orbits are very close 
to the level lines of J^^\ even when /i is of the same order as Pmin- The dynamics on 
the quasi-invariant curve J(s, 6) = Jq is approximated by ipi = 4'o + /^c( Jo)r2( Jq, fi). 



In section |6.2| , we will say a word on adiabatic invariants for billiards in non- 
convex boundaries. Quasi-invariants can also be used to study billiards with piece- 
wise smooth boundaries. As long as the particle hits the same smooth piece of the 
boundary, the corresponding quasi-invariant changes slowly, so that one can estimate 
the outcoming velocity. However, as soon as the particle passes from one smooth 
piece to another one, a jump of the value of J occurs. This is particularly clear in 
the stadium, where J = ^ is even an exact invariant on each straight line and arc of 
circle. The jumps of 6 occurring at the junction points allow the particle to explore a 
large fraction of phase space (see also [p!o[] ). Thus, billiards in insufficiently smooth, 
convex domains do not become integrable in the strong field limit (as suggested in 



6 More chaotic billiards 
6.1 The billiard in a square 

Billiards in convex polygonal domains, in the presence of a magnetic field, can 
be expected to show a chaotic tendency, because we will always be in the regime 

= Pmin < P < Pmax = OO. 

We have concentrated our study to the case of a square of side length 1. This 
billiard is integrable in zero field, where every periodic orbit is parabolic. In fact, 

^Namely, K = /( J^^)) + ©(^x^), where /(J) = - cos J(l + 2 sin^ J) sin^^ J. 



Figure 7: Periodic trajectories in a square. Trajectories of period 4: (a) type I, 
(b) type II, (c) type III, (d) asymmetric trajectory at /i = 0.502, bifurcating with 
type I at yU = |. (e) Trajectory of period 12 which bifurcates with type I when 

H = (f) Stable symmetric trajectory of period Ak, = 2 for = 



the periodic orbits occur in families, which can be indexed by the slope of the 
trajectory, and whose members can be indexed by the arclength of any collision 
with the boundary. 

In low magnetic field, only isolated orbits subsist, some of which are hyperbolic, 
and some elliptic (fig.pa)0. When the field increases, most of the latter vanish or get 
unstable, and structure of phase space is dominated by symmetric period-4 orbits 
(fig.|^-c): type I exists for /x > type II for /i > |, and type III for < /i < |. 

Using (1^), we find that type II and III orbits are hyperbolic, except for /x = 
where type III and type I undergo saddle-node bifurcation. For type I orbits, we 
find that the stability matrix satisfies 



-TiS, = 8t'-8t^ + l, t = -Trgi= ^ \ (51) 

2 ^ ' 2 ' 1 + V8/i^-l' ^ ^ 

so that these are elliptic, except for /i = (saddle-node bifurcation), /i = ^ (bi- 
furcation with 4 elliptic and 4 hyperbolic asymmetric orbits of period 4, fig.^, §3), 
and u = — 7^ — (no bifurcation). 

2V3±2v^ 

An interesting resonance phenomenon occurs when fi = 6*4 is then a cubic 
root of 1, and type I orbit becomes unstable because of a "squeeze effect" bifur- 

^Thc first correction to the zero-field generating function goes like the magnetic flux through 
the trajectory, i.e. the signed area enclosed by the trajectory. For orbits whose period is a multiple 
of 4, it can be seen that this area varies quadratically with the arclength indexing the orbit. 




Figure 8: Phase portraits of the bilhard in a square, (a) /i = 30, (b) /i = 0.502, 
(c) /i = |. "Squeeze effect" bifurcation around = ^ = 1.118. . .: (d) /i = 1.25, 
(e) /i = 1.118, (f) fi = l. 



cation ( ||A1|| , p. 392) with a hyperbohc orbit of period 12 (fig.^e, |^d-f)|. Numerical 
simulations show no other stable periodic orbits in this case, leading us to the con- 
jecture that the billiard in a square is ergodic for this particular value of fi. Figure 
^ illustrates this phenomenon by showing a numerical estimate of the size of regular 
components of phase space, which vanishes for /i = ^ ~ 1.118 . . .. 

There are no orbits of period 4 when < However, stable orbits of higher 
period exist for arbitrarily high values of the magnetic field. Figure |^ shows a 
symmetric trajectory of period 8, doing 2 bounces on each side of the square. We 
can construct similar trajectories of period Ak, k > 1, reflected k times on each 
side. Detailed calculations (see appendix 0) show that such orbits are stable for 
almost all /i in some interval containing but these intervals do not overlap when 
k > 2, leading us to the assumption that mixed and ergodic behaviour alternate for 
decreasing /x. 



6.2 Boundaries with negative curvature 

In Sinai billiards, whose boundaries have negative curvature, it is known that all 
periodic orbits are unstable in zero field (as is suggested by the fact that Tt{DT) < 
—2 if we put X = in (0)). This is no longer true in non-zero field, because the 
curved trajectories, although locally dispersed, may converge again. 

^The same kind of bifurcation affects the small diametral orbit in the ellipse when A = 2, 

VyLt > 1. 



Proportion of regular component 
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Figure 9: Numerically estimated proportion of regular component in phase space 
for the billiard in a square, as a function of fi. Actually, we show the relative area 
of phase space not occupied by one chaotic orbit, with a precision of 4 ■ 10~^. 



Let us first consider billiards outside smooth convex boundaries with extremal 
radii of curvature pmin and Pmax- If we fix a trajectory of the outside billiard, and 
complete every arc of trajectory to a full circle, we obtain its "dual trajectory" as 
the set of all complementary arcs. If the dual trajectory never crosses the boundary, 
then it corresponds to a real trajectory of the inside billiard, since the law of specular 
reflection is satisfied by construction. A sufficient condition for this to be true is 
that any circle of radius /i can intersect the boundary at most twice. In this section, 
let us call this the "/i-intersection property". This property is actually granted in 
two cases: in strong field /i < pmin (as we prove in lemma 3 of appendix 0), and 
in weak field fi > pmax (as can be proven in a similar way). In these cases, to any 
orbit (. . . (s_i, M-i), (so, uq), (si, ui) . . .) of the outside billiard, corresponds the orbit 
(. . . (— si, —Ml), (—So, —Uq), (— s_i, — M-i) . . .) of the inside bilhard, and reciprocally. 
Thus, the billiards inside and outside the boundary are perfectly equivalent, so that 
we can apply our results of section |^ on existence of invariant curves and adiabatic 
invariants (theorem 1, proposition 3 and the corollary of theorem 2). In particular. 



Figure 10: The billiard outside two circles: (a) chaotic trajectory, (b) one of the few 
stable trajectories. 

the billiard outside a circle is integrable. 

For intermediate values of the magnetic field (pmin < /i < Pmax); the dual of 
any outside trajectory will not necessarily be an inside trajectory. However, the 
equivalence is true for some special trajectories, namely those which are sufficiently 
close to M = +1 for the inside billiard (and remain so because of theorem 1). 

Consequently, we predict that if we compare the phase portraits of the billiards 
inside and outside a given smooth convex curve, they will be the same (up to an 
inversion) for /i < pmin and /i > Pmax- When pmin < p < Pmax, the region m ~ 1 of 
the interior portrait will be the same as the region u ~ — 1 of the exterior portrait, 
but other parts of them will be different. 

Let us consider next the case where the boundary dQ has both convex and 
concave parts, but a bounded curvature: |k(s)| < = y—is. We can argue by the 
means of geometrical properties that such a billiard must possess invariant tori and 
adiabatic invariants in sufficiently strong magnetic field, if dQ is smooth enough. It 
should be clear that the interior of the domain has to be connected, since otherwise 
we would have several independent billiards. We can thus define r* such that the 
//-intersection property is satisfied for /i < r* This radius has to be positive but 
can be arbitrarily small (consider billiards shaped like sand-glasses or peanuts). 

In section [5.2| , we already mentioned that the expression (^) of the bouncing 
map near 6 = n can still be valid if /i < Pmin(l — namely if 6 is sufficiently 
close to 71 that the trajectory remains within a distance 2r* of the boundary. Thus, 
invariant curves still exist near 6 = n. Furthermore, we note that the only point in 
the proof of proposition 3 where we use the convexity of the boundary is to show 

^°Note that the bouncing map is continuous iff the /i-intersection property is satisfied. Indeed, 
discontinuities by tangency exist iff one can construct a Larmor circle tangent to the boundary, 
and intersecting it at two other points 



the /i-intersection property. This is actually satisfied for /i < r*, so that we conclude 
that proposition 3 and its consequences on adiabatic invariants remain true in the 
present case, for sufficiently small fi. 

Finally, we also studied the billiard outside 2 circles of radius 1, centered at 
(0, ±A) (fig.|TOD. If A > 1, it is important to note that Mi defined in section |^ 
is not the only integrable component of phase space. Indeed, trajectories with 
u > u* = min|l,2A — 1 — ^^^^-^| always touch the same circle, so that u is a 
constant of motion. In particular, when fi < X — 1, the billiard is integrable. When 
/X > A — 1, the map is discontinuous whenever trajectories become tangent to the 
boundary, and numerical simulations show that the component of phase space u < u* 
is almost filled with a stochastic sea with positive Liapunov exponents. However, one 
can find stable periodic orbits for particular values of the parameters (see appendix 

so that the billiard is not always ergodic. 



7 Proof of the adiabatic theorem for maps 



7.1 Construction of the change of variables 

The change of variables is constructed as the composition of several elementary 
substitutions, that each increase by one the order of the terms containing the phase. 
The number of these substitutions is k if the map is and of order [1 / /i] if the 
map is analytic. 

We assume by induction that after steps, the map reads 

ifQ + /ic(/o) ^n{Iq, /i) + /i^^^aAr(v9o, -^0, At) 

/o + /ic(/o)^ /ieAf(/o,/i) +/i^+^/5Ar(v3o,/o,/i) , 



h 



(52) 



where we write (v^, /) = (v?*-^-*, Z^^-*) for brevity, and VLm{I^iA = ^^o + /^^Af(-^;/^)- We 
can assume that < >=< [3^ >= 0, where < F >= Jq Fd(p denotes the average 
of F. 

The change of variables (</?, /) (tp, J) = {ip^'^^^\ /(^+i)) jg defined by 

^ = y. + /i^+V(¥',/,/i) 
J = I + ^''+'c{I)g{^,I,fi), 



(53) 



where 



•fi 1 

9(^,1,1^)- < 9> (^,yu), 
/ c (/)^o(¥'', ^, /i) - t^(^n{^ , /i) 



(54) 



We will use the following properties: 
1. If F{ip) e then 



F{Lp)dip = /ic(/o) l^Nilo, At)^(v5o) + IJ'ri{(fo, Jo, //)] 
= /ic(/o) [noF{ipo) + /ir2(v2o, /o, y")] • 



(55) 



Indeed, carrying out the change of variables ip = ipQ + ^c{Io)x, one obtains 



V>1 



F{ip)dip = /ic(/o) 



/ F{ipo)dx + fic{Io) / F'xdx+ / Fdx 

Jo Jo iOjv 



2. If F{I) e C\ then 



F(/o) = F( Jo) + /i^+^c(/o)r3(<^o, /o, 



(56) 



where r3 = —g{ipQ,lQ,fi)F'{I). In particular, c(/o) = c{Jq)Ro{ipo, Iq, fi) , where 



Substituting ( |52D in (|^), we get 



Ji — Jo 



Af+l 



'cih)g{^i, /i, /i) - c(/o)^(</?o, /o, /i)] + 



+/i^c(/o)^ [07v(/o, /i) + ^^pNi^o, Iq, Ai) 
Now, using the above properties, we obtain 
cili)9ivi, A, yu) - c{Io)g{(pi, Jq, /i) : 



(57) 



{h-Io)dj[c{I)g{^i,I,fi)] 

/i^c(/o)^i?i(V2o,/o,/^), 
1 /■'^i 

-ficilofpNi^oJo,!^) + 
+/^^c(/o)^i?2(v5o,^,/^), 

/i^+^c(/o)/?3(</'0,/o,/i), 

/x^+^c(/o)r4(v9o,/o,/x). 
Replacing this in (pT]), we finally get 

Ji = Jo + /i^c(Jo)2e^(Jo,/i) + /x^+MJo)^/37v+i(<^o,/o,/i), 



c(/o)^(¥'i, ^, /^) - c(/o)^(¥?o, ^0, /^) 



c(/o)6iv(/o,/i) - c(Jo)6iv(Jo,/i) 

c(/o) - C( Jq) 



(5^ 



(59) 



where /Sat+i = Rl{Ri + R2 + R3) + RqRa and -R4 = 6Arr4. The final step is to write 
PN+ii'^o, lo, f^) = PN+iii^o, Jo, 1^), which can be done for small enough /i, as we shall 
check later, and Qn+i = ©iv + fJ^^~^^ < $n+i >,Pn+i = Pn+i— < Pn+i >■ 
Proceeding in a similar way for ip, we obtain 



-01 - Ipo 



/ic(/o)fi^(/o,/i) l + /i^+ic'(/o)(?o(<^o,/o,/i) + fi'^-^'MW 



,N+3, 



c{Io)flN{Io, IJ) = c{Jo)^n{Jo, /i) + + c(/o) [-flNih, fJ')c{Io)g{(po, Iq, jj) + /i-Re] 

c(/o)c'(/o)fiAr(/o,/i) <g> (/o,/i) = c(Jo)c (Jo)fiAr(Jo,/i) < fi- > ( Jo, /") - IJ,c{Io)R7, 

so that 

t/)! = V'o + hc{Jo)CIn{Jo, h) + /i^+^c( Jo)aiv+i(v9o, Jq, /i). (60) 

where oat+i = RoiR^ + Rq + Rj) and fi^v = f^Af(l — fi^~^^c' < g >). Finally, we 
proceed in the same way as for /5 in order to define f2iv+i, aN+i- 



7.2 Bounds and domains of analyticity 

In the case where the map is analytic, we define 



T{D) = {(v?, /) I \lmip\ , |Im/| < D,-D <ReI <1 + D} 



(61) 



We assume that for /) G F 
bounds |a7v| , \Pn\ < M^, 



N 



Dn) and < /io, (|52D is analytic and satisfies the 
\Qn\ 



< Wn. We introduce numbers D'^ < D'j^ < 
Dn, such that (v?o,-^o) G r(D^) implies {(fi, h), {'ipo, Jo) E T{D'j^). Thus, D'^ has 
to satisfy the condition 



D'^<D'j^~ ko max {/io + /i^VTjv + fio^^M, 



(62) 



From now on, the numbers ki,Xi,Ci will designate constants which are uniform 
in N and /i. Using Cauchy's inequality to bound the derivatives appearing in the 
expressions of the Ri, we obtain that for {ipo, Jq) G F(D^), \Ri\ < Ki, where 



K2 
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(63) 



and the new bounds have to satisfy 

[ i=l i=5 J 

Wn+1 > WN + C2fi^^\WNMN + MN+l). 



(64) 



Finally, we have to consider the effect of the change of variables on the domain 
of analyticity. 

Lemma 1 If F{{p, I, fi) is analytic for (y?, /) G F(D), then F{ip^ J, /i) = F{Lp, /, yu) 



N 



IS analytic for {^p, J) G T{D - Acofio'M, 

Proof: Write x = {ip,I),y = {ip, J). From (|53D , we have y = x + G{x), with 
\Gi{x)\ < Co/i^+^M^. If ^{x,y) = x + G{x) - y and \diGj\ < e, then det(9^.<l>) > 
l-2£-2e'^. liye F(L)-4co/i^+^M^), then x G F(L) - 3co/x^+^Mjv), and Cauchy's 



inequality implies e < 1/3, so that det(c}i.$) > 0. The implicit function theorem 
implies that we can write x = ^{y) where ^ is analytic, and so is F{y, fi) = F{^{y), fi). 



□ 



Hence, we must take 



N- 



(65) 



7.3 Evolution of the bounds with 

We now choose domains of the form = Dq — NdfiQ, D'^ = Dn — d'fiQ, d' < d. 
We wish to show by induction that for a suitable choice of d and d' and /io small 
enough. 



Mn 



(2/io; 



(66) 
(67) 



Replacing this in 

Ki < kiMN/d'no, i = 1 



, we get i^o < [1 - yUoAoMo] ^ < 2 if /iq < [2AoMo] \ and 
. 7. Thus, we can take 



Mn+1 = Cim— 



M, 



N 



d'no 2/io 



(6^ 



where d' = 2cim and m = max |4 1]f=i ki + 2k4, 2 I]I=5 ^j}- 

To prove {^), we take c3 = Wo + C2M0, W = Wo + C3 and /iq < Wq \2c2M0W 
so that 



C2M0 

2N 



< 



2N 



Wo C2M0 
2 2 



C3 



27V+1- 



Finally, it is easy to check that the conditions ( |62D and (]65|) are satisfied if we 
take d = d' + 4coMo + fco max {l + fi*(W + Mq), Mq}. 
The last step in order to prove (^61) is to take 



N = iV(/xo) 



dfxo 



(69) 



In this way, Wn{^io) is bounded, -DAr(^o) still positive so that r(M7v(^o)) 
empty, and the terms containing the phase are bounded by 



(70) 



Finally, we also see that the terms defining the successive changes of variables, 
i.e. Ho^^/n, where ~ M^, add up to a geometrical series converging towards a 
term of order fj,Q. 



7.4 Conservative case 



It is easy to check that if (^) preserves the measure c{I)p{ip, I, fi)dipdl, p{ip, I, 0) = 
1, then (|5^) preserves the measure c{I)pn{'^, I, fJ')d(pdI, pm{(P,I,0) = 1. As long as 
(^) is and c G C^, we can apply the 

Lemma 2 // ( fj^ j preserves the measure c{l)pdipdl, then 0jv(/o,yu) = 0{p^^^) 
andd^p = 0{p^+^). _ 

Proof: We assume by induction that Qiyc{I) = p^Qn and p{(p, I, p) = po{I, p) + 
/i^+^Pi(v9, /, p), which is clearly true for p = 0. We have 



c(/i)p(v?i,/i,/x) 
c(/o)p(v2o,/o,/i) 



det 



1 + p^+2^(/o)9^«Ar 0{p) 
y 0(p^+2) 1 ^ pf+25,(ce;v) + P^+29,(c2/3^ 

1 + pP-^^di [c(/o)ejv(/o, 0)] + (cS^aiv + 5/ \c'p_ 



?7V 



+ 



(71) 



On the other hand, making expansions, we obtain 
c(/i)p(y?i,/i,p) , , .p+2 



1 + c(/o)fio5^Pi(<^o, /o, 0) + c'(/o)e^(/o, 0) 



+ 



+/x^+2c(/o)c'(/o)/3,v(^o, /o, 0) + 0{p^+') + 0(p^+3). (72) 



Comparing (^) and (|72D, we obtain for p < 

^od^pi{(poJo,0) = 9/eAr(/o,0). 



(73) 



The left-hand side is a periodic function of with average 0. Hence, the right- 
hand side must vanish. This implies d^pi = and Qn = const. But this constant 
must be zero since Bjv is divisible by c(/). This proves p = by induction. For 
p = N, one obtains again that OjOn is a periodic function of (f with average 0, 
which proves the lemma. □ 

In the C*^ case, we first do — 1 changes of variables, so that the map takes the 
form (p2D with N = k — 1 and is C^. Then we can apply the above lemma, and 
finally we do the change of variables (|53|) with N = k, which proves (|49|) . 

In the analytic case, after N{pq) changes of variables, the lemma implies that 
we can write Qn{,I,p) = ^J'^^N{I, p), and Cauchy's inequality gives 



OAr(/,p) 



< 



W 



(/^o-l/^l)^' 



(74) 



For pi = po/3, we have 



N 



< 



w 

2^' 



Pi Pn 



< 



Mo 



which proves 



(75) 
□ 



8 Conclusion 



In this work on classical billiards in plane domains, we found that some properties 
of bouncing maps, which are known in the Euclidean case, can be generalized to 
the situation where a magnetic field is applied perpendicularly to the plane. Exact 
expressions for the Jacobian matrix and a generating function help to find peri- 
odic orbits, analyze their stability and compute bifurcation values. Perturbative 
calculations are improved. 

Some aspects of the behaviour of billiards in convex domains are well understood: 
if the boundary is sufficiently smooth, existence of whispering gallery modes prevents 
ergodicity, and when the magnetic field B goes to infinity, the bouncing map behaves 
like a perturbed integrable map. We were able to construct quasi-invariant s, which 
are conserved for a time of order B^~^ if the boundary is C^, and of order when 
the boundary is analytic. Some of these properties remain valid for a more general 
class of billiards, the boundary of which is not convex, but has a bounded curvature. 

Chaotic behaviour seems to be created by two different mechanisms: non-linearity 
(responsible for instance of chaotic components near separatrices) and singularities 
(due to non-smooth boundaries or discontinuities by tangency). Non-linearity alone 
is not sufficient to make a magnetic billiard ergodic, since we showed that there 
always exist invariant curves if the boundary is and tangencies are impossible. 
This implies in particular that billiards with boundaries of negative curvature are 
not necessarily very chaotic. In any case where the billiard is of the mixed type, 
the really challenging open problem remains to prove the existence of a chaotic 
component of positive measure. 

Finally, we discussed two examples of billiards fulfilling necessary conditions for 
ergodicity (i.e. their bouncing maps have singularities). Our study of symmetric 
periodic orbits in the square lead us to the conjecture that ergodic and mixed dy- 
namics alternate when the magnetic field increases. However, only the first value of 
supposed ergodicity, fx = lies in the region > |, where no trivially integrable 
component of phase space exists. The scattering billiard outside two circles shows 
strongly chaotic dynamics, but still possesses elliptic orbits for some values of the 
parameters. 
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A Proof of proposition 1 

We consider the arc of trajectory in figure If 



a = arg [(X(si) - X{so)) + i {Y{s,) - Y{so))] (76) 



is the angle between the x-axis and PoPi, then we have 

00 = a-To-x, 

6*1 = Ti-a-x- (77) 
Now, using (I) and f = (X(si) - X{so)f + {Y{si) - F(so))^ we find 

^ = l,[{Y{s,)-Y{so))X'{so)~{X{s^)-X{so))Y'{so)] 
sin a cos To — cos a sin To sin(6'o + x) 



(7J 



dso 

dx 1 1 di cos(6'o + x) cos(6'o + x) sin X 

dsQ cosx2/i(9so 2/icosx ^cosx ' 

where the last equality comes from the zero-field generating property of £: 

^ = ][-X'{so){X{s,)-X{s,))-Y'{so){Y{s,)-Y{so))] 

= — cos a cos To — sin a sin tq = — cos(6'o + x) • (79) 
Collecting terms, we get 

dOo sin(go + 2x) 06, smOp 

^ ff '^0' ^ "'ff ■ v^^) 

OSq tcosx OSq icosx 

Similarly, we find 

dOo sine, 89, sin(^i + 2x) , 

osi ccosx uSi tcosx 

These quantities give dO^ and dO, in function of ds^ and ds,. Solving a linear 
system, we can express ds, and dO, as functions of ds^ and dO^. Finally, using 
du = sin 9d9, we obtain the equations (0). □ 

B Proof of proposition 2 

Taking the derivative of (O) and using (0), we get 



di V 4/i 
From simple geometry, we obtain 

dA 1 



1- — = cosx. (82) 



2-m(^o + x), (83) 



so that the derivative of (|T3|) is 



dG 1 OA OB. di ^ ^ / .^ 

^ = \- -^Tr~ = -smxsm^e'o + x) - cosxcos^e'o + X) = Uq, (84) 

OSq fJ^OSo Oi OSq 



where we have used ([79|) and (0) again. We proceed in a similar way for s,. □ 







-1 




Figure 11: Behaviour near m = — 1 of the bilhard in an elhpse when pmin < /i < p^ax 
(/i = 1, A = 2): (a) chaotic orbit near u = —1, (b) the preimages Vn, n = 1 ... 15 of 
the hne u = —1. 



C Discontinuities of the map due to tangencies 

The bouncing map is discontinuous either in the corners of the boundary (s G Ei), 
or if the trajectory becomes tangent to the boundary. As can be seen in fig.|b, 
such a tangency means that there are two arbitrarily close initial conditions -21,^2 
such that Tz2 is far away from Tzi but close to T'^Zi. Then there exists a point z* 
between zi and Z2 such that Tz* = (s*, —1) (or (s*, +1) if dQ is not convex). 
Let C = [0, \dQ\) X {—1, +1}. It can be written as the union of 3 disjoint sets: 

• the set Co of z G C such that T~^z does not exist; 

• the set Ci of z E C such that T^^z G C; 

• the set C2 of z E C such that T^^z ^ C. 

The set V of lines of discontinuity of T is T^^C2. 

To construct the sets Ci, we first look for circles of radius fi which are inscribed 
in Q, i.e. circles contained in Q and tangent to dQ at two points or more. Let E 
be the set of abscissas Si of these contact points. Let E be set of abscissas Sj such 
that |p(sj)| = p. Then Co, Ci and C2 are delimited by points with abscissa in E, E 
and E2 (see |[RB|| for a more geometric interpretation of these sets). 



We illustrate this construction in the case of an elliptic boundary (ITsD , using ip 
instead of s. Here we take C = [0, \dQ\) x { — 1}. Inscribed circles exist if p^am = 
A^^ < p < 1, and their contact points are solutions of 



sin^^. = ^^^. (85) 



If Pmin = A ^ < /i < Pmax = A^, there are points y^j, given by 

If they exist, the solutions can be ordered < < < <^2 < <^2 < "^3 < <^3 < 
< ^^4 < 27r. They satisfy 

T-i(<^„-l) = ((^5-^,-1) for z= 1,2,3,4, (87) 
hm r-i(¥;,-l) = (<^„-l) for J = 1,3. (88) 

Table |I| shows the resulting subsets of C. Following ideas in ||DR|| , one can look 
at the successive preimages Vn = T^"C2 of the lines of discontinuity. Numerical 
simulations show that they densely fill the stochastic layer near u = —1 (see fig.pJ]). 

Table 1: The sets Co, Ci and C2 for an elliptic boundary, 
/i G Co Ci C2 

(0,Pmin] C 

(Pmin,!] (^2,^3] U (^4,^1] [^1,^2] U [(^3,^4] {^l,(pi)U {ip2,^2p 

X{-1} X{-1} (^3,^3) U (^4,<^4] X {-1} 

(l,Pmax) (^2,^3] U (^4,^1] (^1,^2] U ((^3,^4] 

X{-1} X{-1} 

[/Omax,Oo) C 



D Proof of proposition 3 

The following geometrical lemma describes some properties of convex plane curves 
(which seem quite obvious if one draws a picture). 

Lemma 3 Let D be a circular segment of angle 2a G (0, tt] and radius 1. Let 
C be a strictly convex curve with extremities on the vertices Aq,Bq of D, and 
making acute or right angles with the chord AqBq (see fig\T^). Then, if the curvature 
of C is everywhere less than 1, it is entirely contained in D and shorter than 2a. 



Proof: We introduce coordinate axes as in figure 12. For — sin a < a; < sin a, C 



can be described by a function y{x), and the curvature is 

This equation can be integrated between the abscissa Xq of the maximum Sq and x: 

k{x) 



y\x) 



2 



k{x) = I K{xi)dxi, |A;(a;)| < |x — Xol . (90) 

Jxo 
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Figure 12: Geometry of the circular segment of lemma 3. 

For xq — l<a;<xo + l, this implies 

\y{x) - (y(xo) - 1)]' + (x - xo)' > 1, (91) 

and hence C lies above the circle of radius 1 tangent to it at its summit 5*0 . Now 
assume by contradiction that C lies outside D between two points P and P' . Let S 
be the point of C furthest away from PP' . The circle of unit radius tangent to C at 
S cannot intersect D twice between P and P', which contradicts (pi]). 

If P and P' are two points of C, we denote by PP' the length of C between these 
points. If P{x) = (x, y{x)), then for xq — 1 < x < Xq + 1, 

S'oP(x) = / a/ 1 + ?/'(xi)2(ixi < Arcsin |x — Xol . (92) 

Jxo 

Assume that Xq < 0. Then AqAi < Arcsin(xo + sin a) + Arcsin(— xq) < « by 
convexity of Arcsin. Considering the circular segment of vertices Ai and Bq, one 
obtains that AiBq can be divided into two parts, one of which is bounded by |. 
Repeating this procedure, we see that C can be divided into pieces whose lengths 
are bounded by a geometrical series, converging towards 2a. D 

Corollary Let C be a plane convex curve, whose curvature satisfies < 
< k(s) < — !— < oo. Let T be a circle of radius u < Pmin- Then 

Pmax ^ ^ Pmin 

1. F can be tangent to C at one point at most. 

2. F can intersect C at two points at most. 

Assume that F intersects C at P and P' and call Co = C (1 IntF. 

3. Co is shorter than TCfi. 

4. The angles between PP' and Co are acute. 

Proof: Assuming that 1. or 2. are false, one can construct a counterexample of 
the lemma (one may have to translate F). 3. and 4. are obvious if one considers all 
the circles of radius fi containing P and P'. □ 



Figure 13: Geometry of the trajectories for a convex billiard in strong magnetic 
field. 



We now proceed to the proof of the proposition. 

Choose two points P and P' on dQ. They can be connected by two, one or no 
arc of trajectory, depending on whether their distance is smaller than, equal to, or 
larger than 2/x. In each of the first two cases, we call sq and si their abscissas, as in 



figure 0. We define 

^ ^ (sq-Si) (mod \dQ\) 

The corollary and our sign conventions imply that < o" < vr. Moreover, the arcs 
cannot intersect dQ at a point different from P or P'. 

Following the proof of proposition 2, it is easy to show that the two trajectories 
can be described by the generating functions 

G^iso, s,) = ~A{s„ So) ± fib , (94) 



in the sense that dG^ = UQdso — ufdsi. Here, £ G (0, 2/i] and A are defined in fig.[T3|, 
and 



b{x) = Arcsinx + xvl — x^ 

b'{x) = 2Vl -a;2. (95) 

The functions A and i are directly related to the shape of dQ, that we describe 
using the function r(s), which has the properties 

1. T{s+\dQ\) = r(s) +277, 

2. t'(s) = k(s) G [^,^1, 



r\dQ\ . , , 
Jo 



If -^(so, Si) is the vector connecting the points of abscissas sq and si, we have 
£{so, Si) = / t{s)ds 

J so 

i{so, si)2 = r ds r ds' (t(s) |t(s')) = r ds r ds' cos [r(s) - r(s')] 

J SQ J SQ J SQ J SQ 

A{so,si) = r dsl\i{so,s) A t{s)\ = l r ds r ds' sin [r{s)-T{s')]. (96) 

J sn Z Z J sn J sn 



:h£{so,s)Atis)\ = l r 

Carrying out the change of variables s = So+fJ.crt, we obtain £(so, Si)^ = /i^cr^/(so, fJ^cr), 
^(■5i,so) = i/i^(T2J'(so,/icr), where 



1-0 1-0 

/(so,/iO") = dt dt' cos T{so^ fiat, ^at') 

rO rt 

/'(so,/icr) = J dt J dt' sin T{so, ficrt, fiat') 

T(so, fiat, fiat') = r(so + fiat) — r(so + fiat'). (97) 

If the boundary is C^, one gets T = 0{fia), 1 = 1 + 0{fia), I' = 0{fia), and 
thus G^{so, si, fi) = fig^{a. So, fi), where g^{a,so,0) = ±^ (f ) does not depend 
on So- From this fact, we can aheady guess the structure of the map. Indeed, 



d^g^ + fids^g^, uf = d„g^, and thus = Mq + fidg^g^ = + ^(/i^). 



However, G is not sufficiently smooth around £ = 2/i to be expanded. We proceed 
in a slightly different way: taking derivatives of (p6D, we obtain 



where 



^^ fi dA fi OA fi 
aso t osq 2 osi 2 



J{so,fia) = J dt cos T{s 0,0, fiat) 
K{so,fia) = / dt sin T{s 0,0, fiat) 

rO 

K'{so,fia) = J dt cos T{so, fiat, —fia). (99) 

Note that if the boundary is C'', the above integrals are all C^^^. Differentiating 
(0), we get 



, dC^ & di IdA 1 /i 1 



1 A 1 J' 



UO + -.K) = ^l--a^I^ ^ 

AI{l-ul) =AIsin^eo = a'^{K'^ + J^) I - Aa cos 9 oK I + A{I - J^) 

= $o(a,/i,So,^^o). (100) 



When the boundary is C^, we have 1=1 + 0{jj,'^a'^), J = 1 + O^^i^a"^) and 
K = ^K,{so)]icr + 0(/i^cr^), so that $o = (x^ [1 — 2 cos6oK,{so)fi + 0{fi^)]. 
If we write a = 2 sinOoi], then 



= ^{r],fi,So,9o 




1—7] 1 — COS 6'oK(so)yU + 



(T=2 sin^o'? 



(101) 

The function $ has the properties ^{7]q,0, so,9q) = 0, where 7]q = 1, and 
dri^{7](), 0, So, 6*0) = —1. Thus, the imphcit function theorem imphes that in a neigh- 
borhood of /i = 0, ?7 can be expressed as a C^~^ function of /i, sq and Oq. This 
function can be constructed using Newton's method: if ?7o = 1 and ?7„+i(so, 6*0, /i) = 
r7„(so, 6*0, /x) + $(?7n(so, 6^0, /i),/i, 30,6*0), then $(?]„, /x, so, 6*0) = 0(/i"+^). For /c = 3, 
we get 

7]{so, 00, yU) = 1 + cos eoK{so)]i + 0{]I^). (102) 

The first equation of (|l|) is obtained by using si = sq — ficr = sq — 2/i sin 6*077, i.e. 



a{so,9o,]i) = -{l-7]{so,9o,]i)) eC'' ^. 
]i 



(103) 



Note that if we had used the variable u instead of 6, we would not have been able 
to apply the implicit function theorem when u —^ ±1. Indeed, the map expressed 
in the variables (m, s) contains the factor \Jl — Uq. 

To obtain the second equation of (^Tj), we observe that (see (^)): 



Oi + X)- (^0 + x) 
TT — Arcsm -— — 

\t OSi 

K \ 



7T — Arcsin 



2dA\ 



K' 



Arcsin ( \ — Arcsin ( \ = 0(/i^cr^^ 



(104) 



where we have used point 4 of the corollary to choose the determination of Arcsin. 
By consequence. 



b{so,9o,fj.) 



]i^a'^ 



Arcsin 



— Arcsin 



e C 



k-3 



(105) 



(7=2 sin 9ori 



Higher orders of the expansions of a and b can be obtained by computing expan- 
sions of /, J, K, K' and 7] up to the desired order, and replacing them in ( |103D and 

m. □ 




E Symmetric orbits of period 4/c in the square 

We want to study existence and stability of 4/c-periodic trajectories, similar to the 
one shown in figure where k = 2. As shown in figure |n|, such trajectories can be 
characterized by numbers x, y or x', y', satisfying 



(106) 



( 2kx + 2{k+l)y = l iO<y<l 

\ 2y'^ + 2xy + x^ = fx"^, \ < 2y + x < I, 

j 2kx' + 2{k - l)y' = 1 f < 

\ 2t/'2 + 2xy + x'2 = /x^ jo<a;'<|. 

These equations have 3 different kinds of solutions: 

l+R _ / 1 

^ ~ ^ /i,- < /i < 7^, (107) 



2(A;2 + 1)' ^''-^-2k' 

-l + R 1 1 

— < /i < 



2(fc2 + i)' 2k -"^ V2{k-iy 



where 



1 



i?=^l-(P + l)(l-4A;V), = ___=. 



Using (|^), we find that the stability matrix for k bounces is 



/C ^,{1-C')\(1 -2f,cY-'_(C Ml - (2A: - 1)C2) \ 
' ~ [-i C jl^O 1 ; {2k-l)C j' 

t = ]-TiSk = Ck, (109) 



where 



C = cotg^^ = ^ = -^^. (110) 

y + X y + X 



Now, since 

-^4^ = ^^', ^TTS,k = 8t'-8f + l, (111) 

the total orbit is hyperbolic if \t\ > 1, parabolic if \t\ = 0, 1, and elliptic otherwise. 
Applying this to (|109|) , we find that |/+-orbits are hyperbolic as soon as /i > /i^, 
y_-orbits are never hyperbolic, and ^/'-orbits are hyperbolic if /i > /i^, where 



d = (112) 

We have thus obtained that symmetric 4/c-periodic orbits may be stable only if 
/i^ < /i < /U^. These bounds have the properties 

1 + 

^'^ (113) 
f^k-1 > /^fc if k>2. 

We see that when < /i^, no 4A;-periodic orbit exists. At /i = /i^, a pair of such 
orbits with opposite stability appears in a saddle-node bifurcation. The stable one 
loses stability at fi = fi^. Numerical simulations show that new stable orbits are 
created, but they quickly loose stability for some n = n\. For k > 2, this happens 
long before a 4(/c — l)-periodic orbit appears at /x = fi^_i, and since no other stable 
orbits can be found in the interval, we are lead to the conjecture that the billiard in 
a square is ergodic when fil < fi < 



F Elliptic orbits of period 6 outside 2 circles 



We want to show that the billiard outside 2 circles described in section p.2\ possesses 
elliptic orbits for some values of the parameters. 

The trajectory depicted in figure |TU|b for /i = A = 1 turns out to be linearly 
marginally stable for these values of the parameters if we apply (0). Thus we have 
to analyse its stability for nearby value of and A. 



The trajectory can be characterized by two angles 6' and (p, as in figure |T5|. The 
other angles are then given by Xi = | + 6'' — |, X2 = 6' + ip — ^. Using the relations 
/xsinxi = sin and cos 99 + /xsinx2 = A, we obtain the system 



1 + cos 95 = /i^(l + cos((/) - 26')) 
cos if — fi cos{ip + 9') = A. 



(114) 



For small values of £ = /i — 1 and 5 = A — 1, it has the solution 

e' 



TT 



(115) 



Figure 15: Geometry of period 6 trajectories outside 2 circles. We show a quarter 
of the trajectory, which is symmetric with respect to the coordinate axes. 

where 0(2) stands for terms of order e^, 5^, e5. The Jacobian matrices Mi and 
M2 for the two types of bounces can now be computed using with 6i = tc — 6' 
and iicosxi = /^sin(2xj). The stabihty matrix of the orbit is given by 5*6 = {S^Y, 
S3 = MIM2. We obtain 

t = ^TtS3 = -l-12e + 66 + 0(2), 
^^Tr^6 = 2t2-l = 1 + 485-245 + 0(2). (116) 

Hence, for small S and e, the orbit is elliptic for 6 > 2e + 0(5^), i.e. for A > 
2/x — 1 + 0((/x — 1)^), in a neighbourhood of A = /i = 1. 
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